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^ Abstract 

Consider the set of points formed by the integer n x n grid, and the process that 
in each iteration removes from the point set the vertices of its convex-hull. Here, we 
prove that the number of iterations of this process is 0(n 4 / 3 ); that is, the number of 
convex layers of the n x n grid is Bfn 4 / 3 ). 
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Figure 1: A point set and its decomposition into convex layers. 



1 Introduction 

For many algorithms, the worst case behavior is rarely encountered in practice. This is be- 
cause the worst case behavior might require a degenerate and convoluted input. To address 
this gap between the worst case analysis and a real world behavior, a considerable amount 
of research was spent on analyzing algorithms and discrete geometric structures under cer- 
tain assumptions on the input, including (i) realistic input models [dBKSV02j, (ii) fatness 
[AdBESlT] . (iii) randomness, etc. 



Random points. There is a significant amount of work on the geometric behavior of 
random point sets [R5631 |Ray70[ IWW931 IBar08l lUBSCOOl [JN04] . The question of how the 
Voronoi diagram or the convex-hull of a point set randomly generated inside a convex domain 
behaves had received considerable attention. In particular, it is known that for a set of n 
points chosen uniformly in the unit square, the expected complexity of the convex-hull is 
O(logn), and 0(n 1,/3 ) if the domain is a disk (this bound holds for any convex shape). 



Grid points. Surprisingly, the known results on uniformly sampled points match the re- 
sults known for the grid point set. For example, the number of vertices of the convex hull of 
any subset of the y/n x ^/n grid is O (n 1 / 3 ) , which matches the bound for the random points. 
This phenomena holds for many similar scenarios, see the survey by Barany |Bar08j . 



Convex layers. The decomposition of a point set into convex layers is one possible way 
to measure the depth of a point inside the point set. Formally, the convex depth of a 
point p in a point set P is d p (P) = 1 if p is a vertex of the convex- hull of P, and it is 
d p (P) = 1 + d p (P \ V(CH(P))) otherwise, where CH(P) denotes the convex-hull of P and 
V(CH(P)) denotes the set of its vertices^] This partitions the point set into convex-layers, 
as depicted in Figure [I] In particular, if the points rise out of physical measurements (that 
might contain noise), a point with large convex depth is unlikely to be an outlier. This is one 

1 A point of P is a vertex of the convex- hull only if it is a corner of the convex- hull. Formally, p is a vertex 
of the convex-hull of P is CU{P) ^ CU{P \ {p}). 
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possible definition of robust statistics for points, although this definition has its limitations, 
see |RS04j for details. In particular, Chazelle |Cha85] provided an O(nlogn) time algorithm 
for computing all the convex layers for a set of points in the plane. 

For a set of n points picked uniformly inside a bounded convex domain in IR d , it is known 
that the expected number of convex layers is e(n 2 ^ d+1 ^) |Dal04j . 



Our results. In this paper, we are prove that the number of convex layers of the nxn grid 
is 0(n 4//3 ) . This bound is quite surprising - indeed, as demonstrated by Figure [2 the peeling 
process starts out quite slowly, the first three layers having 4, 8, 8 vertices (independent of 
the value of n), respectively. A priori, it is not clear why this process accelerates and contains 
more vertices. Furthermore, the maximum numbe r of vertices in convex position in an n x n 
grid is 0(n 2/3 ) (th is is well known, see Lemma 2.lh. Namely, somewhat surprisingly, a 



constant fraction of the layers have asymptotically maximum size. Our result matches the 
known result for random points. Note, that although the bounds are similar, the proof for 
the random point set does not carry over to the grid case. 

We also observe that the number of convex layers is Q(n 2 ) if the grid of n x n points 
is allowed to be non-uniform (instead of the integer grid used above). Naturally, in this 
construction, where every point is on two lines where each has n points. 



2 Peeling the grid 

Let Po = G n = {1, . . . , n} 2 , be the nxn integer grid. In the ith iteration, consider the 

convex-hull Cj = C"H(P,;_i), for i = 1, Let Vj, be the set of vertices of Cj. Naturally, we 

consider a grid point to be a vertex only if it is a corner of the convex-hull, and as such grid 
points falling in the middle of edges of Ci, are not in V{. Now, let Pi = P{-\ \Vi. In words, 
we start with the nxn grid, and peel away the vertices of the convex-hull, and we repeat 
this process till all the grid points of G n are removed. Let r(n) be the number of iterations, 
till Pi is an empty set. Here we are interested in the behavior of T{n). See Figure [2] for an 
example of how the generated polygons look like. 

2.1 A lower bound on r(n) 

The following is well known, and we include a proof for the sake of completeness. 

Lemma 2.1. Given any convex set C in the plane, it can have at most 0(n 2//3 ) vertices of 

Gn- 

Proof. Consider a convex set C such that all its vertices are points of G n . The perimeter of 
C is at most An. The number of edges of the convex hull of C of length at least (or equal 
to) fi is at most An/fi. The number of edges having length smaller than \i is bounded by the 
number of integer points of distance at most /i from the origin, and this number is bounded 
by {2ji + l) 2 = 0{ji 2 ). As such, the number of vertices of C is at most O {n / ft + /x 2 ) . Setting 
jj, = [n 1 / 3 ] then implies the claim. □ 
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Figure 2: The polygons generated while peeling the 11x11 integer grid. 
As such, \Vi\ = 0(n 2 / 3 ), which implies immediately that r(n) > n 2 / maxj \Vi\ = f2(n 4//3 ). 
2.2 An upper bound on r(n) 

An integer vector (x,y) is primitive if gcd(x, y) = 1. For an integer fi, let V M be the set of 
all primitive non-zero 

integer vectors (x,y), where < y < x < fi. The following is well known, and we sketch 
a proof for the sake of completeness. 

Lemma 2.2. We have |V M | > c/i 2 ; for some constant c > 0. 

Proof. For a fixed x, consider the vectors (x,y) in V M , such that y < x, and gcd(x,?/) = 1. 
The number of such vectors is the number of integer values of y that are relative prime to x, 
and this number is the Euler's totient function As such, |V M | > Y^i=i 0(0 — ■ the 

last step follows from known bounds, see |HW65] . □ 

In the following, we pick /i to be smaller than n/4, and n is sufficiently large. 

For every vector v G V M , consider the set L v of all lines having direction v that intersect 
the grid points G n . Every line in L v contains at most 1 + [(n — l)/v x \ points of the grid 
(and most lines in this family contain at least |_(?t. — 1)/ points of the grid (the only 
problematic lines are the ones that have short intersection with the square [l,n] 2 because of 
the corners)). 

Claim 2.3. For n > 10, fi < n/4 and v G V M , we have that \L V \ < 4n/i. 

Proof. A line i G L v that intersects G n has an intersection of length at least n with the 
enlarged square [l,2n] 2 . Specifically, the projection of the intersection on the x axis has 
length at least n. Since i has direction v and it contains a grid point, it follows that it 
has grid points on it, that are of distance \\v\\ from each other. On the projection, the 
distance between these points is v x . As such, this intersection contains at least 1 + \ n/v x \ > 
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(A) 



(B) 



(C) 



Figure 3: (A) An active direction v, and the set of lines L v . (B) An inactive iteration for v. 
(C) The next iteration - the two "old" tangent lines no longer intersect the current convex 
layer. 

n/fi points of the grid G^n on it- in particular, the number of such lines can be at most 



Since the lines of L v cover all the grid points of G n , and the vertices of Cj are grid points, 
it follows that L v always contains two lines that are tangent to Cj. If these two tangent lines 
intersect dCi along an non-empty edge, then v is active at iteration i (i.e., v is not active 
if the two tangents touch Cj at a vertex). 

In the following, we slightly abuse notations and use L v fl Cj to denote the set of all lines 
of L v that have non-empty intersection with Cj. 

Claim 2.4. If v is not active at iteration i, then \L V fl C + i| < \L V fl Cj| — 2. 

Proof. If v is not active at iteration % then a tangent i to Cj from L v intersects Cj only at 
a vertex. But this vertex is being removed from the point set when computing Pj+i. In 
particular, the line £ no longer intersects Cj + i. The same argument also applies to the other 
tangent. This is demonstrated in Figure |3j □ 

Claim 2.5. Throughout the process, for a vector v E V M; it can be inactive in at most 2n\i 
iterations. 



An 2 



/{n/fi) = 4n/i. 



□ 
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Proof. Every time v is not active, the number of lines of L v that intersect the active convex- 



hull decreases by two, by Claim 2.4 By Lemma 2.3 there are at most 4n/i lines in the set 



L v , and as such this can happen at most Anfx/2 times. □ 

If the process continues more than M = 4nfi iterations then every vector in is active 
in at least half of the iterations. In particular, if is the number of active directions at 
iteration z, then we have that 



w 

a = 

i=i 



rij > 2n/i |V M | > 2cn/i 



by Lemma 2.2 



Observe, that if n« vectors are active at the zth iteration, then the convex hull of C, has at 
least 2rii edges (and thus vertices) at iteration i. As such, if we set [i = |~n 1//3 / c 1//3 ] = 6 (n 1 / 3 ) , 
we have that the total number of vertices of the convex hulls in the first M iterations is at 
least 

2a > Acnfi 3 > An 2 , 

which is a contradiction, as the initial grid set has at most n 2 points. We conclude that 
the algorithm must terminate after M = An^i = 0(n 4,/3 ) iterations. We thus proved the 
following. 

Theorem 2.6. Starting with the grid G n , consider the process that repeatedly removes the 
convex-hull vertices of the current set of vertices. This process takes 9(n 4//3 ) steps. 



3 Lower bound of Q(n 2 ) for a non-uniform grid 

This section is devoted to describing a set M of n 2 points in the plane where the peeling 
process takes Q(n 2 ) steps. For simplicity assume that n = 2k for some integer k. 

Take a collection of k squares Si, . . . , Sk where Si has length of its side 3 l and the squares 
are positioned such that their centers coincide with the origin. Let L be the set of 4k lines 
that are obtained by extending the segments of the squares into lines. Finally, let M be 
the set of all intersections of lines in L. Notice that each line contains 2k points and that 
\L\ = Ak 2 = n 2 . See Figure H 

Let the peeling process partition M into convex sets Ci, C2, 

Claim 3.1. For every Ci exists Sj such that Ci C Sj. 

Proof. Let j be the largest index such that CidSj 7^ 0. Notice that Ci is centrally symmetric 
as M is centrally symmetric and this property is preserved by the peeling process. If |Cj n 
Sj\ =4 then Ci fl Sj are the four corners of Sj and thus \Ci\ = 4 as Cj is strictly convex. 
Hence |Cj fl Sj\ = 8 and Ci contains points on both vertical and horizontal lines of Sj in 
every quadrant. Let D be the square with corners being intersections the axis and CH(Sj). 
See Figure [5] on the left. Notice that Si C D C CH(Ci fl Sj) for every I < j. Therefore 

d = d n Sj c Sj. ' □ 
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Figure 4: A point set where the peeling process requires Q(n 2 ) steps. 



The previous claim implies that |Cj| < 8 for every i. Hence the peeling process needs at 
least n 2 /8 = Q(n 2 ) steps. 



4 Conclusions 

The most natural question left by our work, is the prove similar bounds in higher dimensions. 
This seems quite challenging, and we leave it as an open problem for further research. 

Let us also note for an interested reader that, according to experiments, the layers in the 
peeling process are getting close to circles as the process is advancing. 
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